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Abstract. Locally periodic rods, which show approximate invariance with respect to translations, 
are constructed by joining N unit cells. The spectrum then shows a band spectrum. We then break 
the local periodicity by including one or more defects in the system. When the defects follow a 
certain definite prescription, an analog of the Wannier-Stark ladders is gotten; when the defects are 
random, an elastic rod showing Anderson localization is obtained. In all cases experimental values 
match the theoretical predictions. 
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INTRODUCTION 

Due to similarities between the stationary Schrodinger equation and the time indepen- 
dent wave equation, several classical analogs of quantum systems have been analyzed 
in the last few years. Among them, mechanical [1], microwave [2], optical [3, 4, 5, 6], 
and elastic systems [7] have been studied along these lines. This approach is sometimes, 
jokingly speaking, referred to as megascopic physics. 

In this paper we shall present a review of some analogs that we have recently discussed 
using one-dimensional elastic rods. Normal-mode amplitudes and frequencies, which 
are the classical analogs of quantum wave functions and eigenenergies, respectively, are 
obtained experimentally and compared with numerical calculations. The experimental 
set up is shown in Fig. 1 where an electromagnetic acoustic transducer (EMAT), con- 
sisting of a permanent magnet and a coil, is used to excite and detect the vibrations of 
the rod. The transducer operates through the interaction of eddy currents in the metal- 
lic sample with a permanent magnetic field and an oscillating one. This EMAT is very 
flexible and selective; it can excite or detect either compressional, torsional, or flexural 
vibrations according to the relative position of the permanent magnet and the coil, as 
shown in Fig. 2. Experimental details are discussed in Refs. [7, 8]. 

We shall first consider the building up of a rod which shows local translational 
invariance. Since we construct it with a finite number N of unit cells (see Fig. 3), we 
refer to it as a locally periodic rod [9]. A band spectrum emerges when N grows. This 




FIGURE 1. Experimental setup to measure the elastic vibrations of an aluminum rod. In the inset the 
EMAT is configured for torsional waves. The exciter should be configurated in the same way. 




FIGURE 2. Different configurations of the permament magnet and of the coil, excite different modes: 
(a) y (b) compressional; (c) y (d) torsional; (e) flexural 



result is also valid for flexural vibrations, which obey instead a fourth order differential 
equation [10, 11]. If the unit cell consists of two small different rods, an elastic analog 
of the diatomic chain results: optical and acoustical bands are formed as N grows. 

We then induce a symmetry breaking by changing the length Iq of one of the central 
small rods that form the locally periodic system of Fig. 3. A normal-mode frequency is 
found in the forbidden band and corresponds to a localized state. The number of rods 
with different lengths is then increased until all of them are distinct. We proceed in 
two different manners. If the lengths of the rods are altered following a systematic law, 
the analog of the Wannier-Stark ladders can be obtained. If the lengths are selected as 
random numbers, the elastic analog of Anderson localization is achieved. 
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FIGURE 3. Locally periodic rod. The unit cell consists of a small cylinder of length Iq and two smaller 
cylinders of length e/2. 

In what follows we shall present our results with the aid of a set of figures, omitting 
the details, which the interested reader can find in the references given to our published 
work. 

BUILDING TRANSLATION SYMMETRY 

In Fig. 4 we present the experimental spectrum of compressional waves of a rod as the 
number N of unit cells is increased. When N is large enough the spectrum consists of 
allowed and forbidden bands, each consisting of N levels. The wave amplitudes are also 
obtained and compared with the ones obtained numerically using the transfer matrix 
method [7]. The waves are extended along the rod. This is consistent with Bloch's 
theorem and what was found numerically for a Kronig-Penney model [12] or for a 
number of other systems [9]. The same result holds for torsional vibrations [7] and, 
even more, for flexural vibrations [11]. This is interesting in itself since the latter obey 
an equation which is not of the Helmholtz type, but rather a fourth order differential 
equation, either the Bernoulli-Euler one or the more accurate one at higher modes, the 
so called Timoshenko equation [10]. The band spectrum for flexural vibration is shown 
in Fig. 5. 

DESTROYING THE SYMMETRY 

We start by introducing in the locally periodic system depicted in Fig. 3 a single small 
rod with a different length I = Iq + h. This is referred to in solid state physics as a 
topological defect, since it destroys the long range order of the locally periodic or 
symmetrical rod. We have measured and calculated the spectrum for several values of 
h [13]. As given in Fig. 6, there appears a frequency in the first forbidden band of the 
ordered rod. This is consistent with what one learns from the band theory of crystals. 
However, in the second forbidden gap two levels lie. In some of the higher gaps two 
frequencies are also found. These state are localized, as seen is Fig. 7. The existence of 
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FIGURE 4. Spectrum of compressional waves of a locally periodic rod. As the number of unit cells 
increases, a band spectrum emerges. 
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FIGURE 5. Bloch theorem applies independently of the structure of the Lagrangian. Band spectra are 
also observed for torsional waves and even for flexural vibrations. Here we show the result for flexural 
vibrations. The circles correspond to the uniform rod, experiment empty circles, Timoshenko's beam 
theory is given by the filled circles. The squares show the results for the locally periodic rod, 12 cells, 
experiment empty squares, Timoshenko's beam theory is given by the filled squares. 



such localization was to be expected from the independent rod model because, if there 
exists a normal-mode frequency proportional to (lo + h)^ 1 , the neighboring rods of the 
defect do not resonate. 
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FIGURE 6. One or more normal-mode frequencies appear in the forbidden gap. 
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FIGURE 7. The states in the forbidden band are localized at the topological defect. Mode 24 continuous 
line, and mode 25, dashed line. 



We can now go on introducing more defects in the locally periodic rod. We first do it in 
an ordered way to simulate, for example, the effect of an external constant static electric 
field on a charged quantum particle that is moving in a periodic potential. This is the 
case, indeed, with the Wannier-Stark ladders (WSL). In 1960, Wannier [14] suggested 
that a static external electric field would destroy the electron band spectrum and instead 
produce a picket fence one, in which the spacing between the energy levels is a constant, 
proportional to the electric field and the lattice constant. The WSL in solid state were 
finally observed decades later using superlattices [15]. They have recently been found in 
optics as well [3, 4]. As will be shown, by modifying appropriately the small rods in a 
systematic way one can obtain an elastic analog of WSL [11]. 

Let us first use an independent rod model. The normal-mode torsional frequencies /j 
of a rod of length /; and wave velocity c, are given by the well known expression [10] 



f. = SL j 



(1) 



where j is the number of nodes in the wave amplitude. We consider two options: 
either / ? - = Iq/ ( 1 + iy) and c,- = a/ G/ p, with y an adimensional parameter, or /,• = Iq 
and C; = c(l +iy). This happens for systems A or B, respectively, which are shown in 
Fig. 8 [16]. Here p is the mass density and G the shear modulus of the rods and c the 
wave velocity. 




FIGURE 8. Elastic rods showing the Wannier-stark ladders 

In the first case, which is the elastic analog of the optical system discussed in Ref. [3] 
we have 

and the difference A/j = — /j is equal to 

n / G 



(2) 



A ^ = 2L\/p J ' r=A/j = J ' Al 



(3) 



which is independent of the index i. 



We should remark that in both cases Aj, is independent of i, a fact which will 
lead to a WSL as we will now show. In Fig. 9 the WSL measured for system A are 
compared with the numerical values obtained using the transfer matrix method. The 
normal-mode amplitudes are localized as one would expect from the independent rod 
model. This is shown in Fig. 10 for several modes in a given ladder; experimental 
wave amplitudes fit very well the theoretical ones, as exemplified in Fig. 11. The same 
conclusions can be obtained by analyzing system B; results for system B can be found 
in References [16, 17]. 
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FIGURE 9. Wannier-Stark ladder spectrum for system A 

A very different effect appears if one now introduces random defects: a disordered 
rod is built. This is reminiscent of the well known Anderson localization. The quantum 
mechanical states in a disordered chain of potentials show always a localization length 

Let us now consider what happens in the elastic case, when we analyze a 1-D system 
formed by N small rods each of length /,-, with = Iq + £ ; and £, an independent random 
variable distributed according to P(Ei) = 1/2 for — 1 < £,• < 1 and P(£i) = otherwise. 
The disorder is an uncorrelated one in this case. 

We first compute the normal-mode amplitudes. It is found indeed that the Ander- 
son transition takes place, the wave functions being localized, as the examples given in 
Fig. 12 and 13 show. Nevertheless, the average localization varies with frequency; see 
Fig. 14. One can also study the statistical properties of the disordered rod spectrum. Us- 
ing, as is common in quantum chaos [18], the nearest-neighbor spacing distribution, we 
find that this distribution also changes with frequency in the elastic case. To characterize 
it we use the Brody repulsion parameter [19]: 



Pp(z)=Az p exp(-<*z P+1 Y> z = s/d 



(4) 
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FIGURE 10. Wannier-Stark ladder wave amplitudes obtained with the transfer matrix method. 
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FIGURE 11. Experimental wave amplitude of the Wannier-Stark ladder and its comparison with the 
theory. 



with A = + a 



13+1 



s the space between frequencies and d the 

mean level spacing. When /3 = a Poisson spectrum is dealt with, for /3 = 1 the Wigner 
distribution, which shows level repulsion, is obtained; the Gaussian orthogonal ensemble 
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FIGURE 12. Localized wave function calculated with the method of the Poincare map of a disordered 
rod. 
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FIGURE 13. Measured localized wave amplitude of a disordered rod. 



(GOE) corresponds to the value /3 = 0.953. 

For the disordered rod the parameter /3 is a function of frequency, showing the rich- 
ness of elastic systems; this is shown in Fig. 15. If one eliminates the frequency, Fig. 16 




io- 41 1 ^ — 1 1 ' 5 u 

5x10 10 2 x 10 

frequency (Hz) 

FIGURE 15. Repulsion parameter as a function of frequency 



is obtained, showing that the repulsion parameter is to a good approximation a linear 
function of the localization length. A similar result was obtained numerically using a 
tight-binding Hamiltonian with nearest neighbor interaction and diagonal disorder [20]. 
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FIGURE 16. Repulsion parameter vs localization length 



CONCLUSIONS 



We have discussed, both from experimental and numerical points of view, different 
configurations of elastic rods. They serve as classical analogs of 1-D quantum systems. 
We first build up a locally periodic rod and a band spectrum emerges. We then destroy the 
local translational symmetry by introducing defects in the rod. If one does it following 
specific rules, the analog of Wannier-Stark ladders can be obtained. If the defects have 
random properties, features of the Anderson transition take place. In all cases, the 
experimental and theoretical results coincide. The elastic analogs are useful to better 
understand several features of one-dimensional models dealt with in the quantum theory 
of solids. In particular, the wave amplitudes can be measured, which is only partially 
possible for microscopic systems. 
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